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The  above  research  papers  can  be  roughly  grouped  into  the  following  areas: 


A.  Applications  of  function  theory, 

B.  The  use  of  summability  methods  and  approximate  conformal 
mapping  techniques  in  the  study  of  iterative  methods. 


We  now  focus  on  area  B,  and  the  associated  paper  by  M.  Eiermann,  W.  Niethammer, 
and  R.S.  Varga,  listed  as  #5  above.  j 

The  Use  of  Summability  Methods  and  Approximate  Conformal  Mapping  Tech- 1  on. 
niques  in  the  Study  of  Iterative  Methods.  — 

To  iterative'y  solve  the  matrix  equation  (in  fixed  point  form)  '  -  - 

.>t  rlb  .tl  on/ 

(1)  x  =  Tx  +  c,  t  Codes 

jA v n 1 1  and/or 


LH.it  1  Tuocial 


□  □ 


(where  1  is  not  an  eigenvalue  of  T),  one  takes  as  the  basic  iteration 

(2)  xffl+1  =  Txm  4- c  (xo  =  a,m>0), 

which  is  well-known  to  be  convergent,  for  arbitrary  a,  iff  the  spectral  radius,  denoted  by 
p(T),  satisfies  p(T)  <  1.  Associated  with  (2)  is  a  semiiterative  method  (SIM),  defined  by 

m 

(3)  ym  =  ^m.i  xi 

«=0 

where  the  (complex)  numbers  {*»«,•  }il’o?m=o  define  an  infinite  triangular  matrix 

*0,0 

u)  p=  T,’°  Xw  ° 

'  ’  ft  2,0  ft?, I  ftl.1 


where  we  assume 

(5) 

On  setting 

(6) 


Y  ^.i  =  1 


?m(z)  :=  Y  7r”‘.<'2’ 


(m  >  0). 


(m  >  0), 


then  pm(z)  is  a  (complex)  polynomial  of  degree  at  most  m,  with  pm(l)  —  1  (m  >  0). 

If  X  denotes  the  unique  solution  of  (1)  (since  1  is  not  an  eigenvalue  of  71),  then  the  error 
vector  em  associated  with  the  iterates  xm  of  (2),  i.e.,  :=  x  —  xm,  satisfies 


(7)  em  =  Tem^  =  ...  =  Tme0 
Analogously,  if  em  :=  x  —  ym,  then  em  satisfies 

(8)  im  =  Pm(T)£0. 


(m  >  0). 


We  now  assume  that  the  set  of  eigenvalues  of  T  (denoted  by  cr(T))  is  contained  in  a 
compact  set  fl  (where  1  ^  0).  The  assumption  of  such  a  set  $7  is  typical  in  applications  of 
iterative  methods  applied  to  matrix  equations  arising  from  physical  problems.  Then,  given 
f 2,  choose  any  infinite  matrix  P  of  (4),  (which  then  defines  the  semiiterative  method  of  (3)). 
Then,  the  asymptotic  convergence  factor ,  with  respect  to  P ,  is  defined  by 


K{U,P)  :=  Jirn  {max  |pm(;)|j 
m—oo  zgn  J 


and  we  set 

(10)  K (fi)  :=  inf  (A  (fl,  P)  :  P  induces  a  SIM  of  the  form  (3)}  . 

Then,  K(Q)  is  called  the  convergence  factor  for  fi,  and  any  P  for  which 


is  called  an  asymptotically  optimal  SIM  (an  AOSIM)  with  respect  to  fi.  The  whole  object  i' 
to  determine  AOSIM 's  for  a  given  compact  set  fi  (containing  the  eigenvalues  of  7  ).  for  such 
SIM’s  give  the  fastest  asymptotic  convergence  rates,  and  these  are  thus  preferred  for  actual 
computation! 

To  describe  our  main  results,  let  C  :=  C(J(oc)  denote  the  extended  plane,  and  let 


(12) 


1M  :  = 


(fi  C  f  :  0  is  compact;  1  6  C\fi;  0°  is  simply 


connected,  and  fi  contains  more  than  one  point}. 

Then,  by  the  Riemann  Mapping  Theorem,  there  is,  for  each  0  €  JM,  a  conformal  mapping 


(13) 


f  *  :  C\{u>  :  M  <  1}  -*  C\(l, 
I  'I'(oo)  =  oc  :  t|>'(oo)  =:  j(fl) 


and 
>  0. 


Now,  there  is  a  unique  u>  with  ty(i2>)  =  1,  and  if  f)  :=  |  -1  >  1  ,  then  this  number  r)  is  the  key 
to  finding  an  AOSIM  with  respect  to  fi  6  IM. 

More  precisely,  we  have  (cf.  [5,  Corollary  [2])  the 


Theorem.  Given  fi  6  IM,  then  any  SIM  generated  by  a  P  satisfying  (5)  is  an  AOSIM  with 
respect  to  fi  iff 

(14)  A'(Q,P)  =  i. 

Thus,  for  any  0  €  IM,  an  AOSIM  always  exists  for  ST,  and  such  a  SIM  can  be  obtained 
from  the  conformal  mapping  of  (13).  What  is  very  interesting,  in  our  opinion,  is  that  there 
is  a  strong  connection  with  the  theory  of  Faber  polynomials  and  the  construction  of  nearly 
optimal  semiiterative  methods,  for  each  fi  6  IM.  To  be  precise,  let  us  normalize  to  *)( fi)  =  1 
in  (13),  so  that 

OO 

(15)  ^(u;)  =  uj  +  ^  oc*  uj~k  ( ju;|  >  1 ). 

k= 0 

If  <f>(z)  denotes  the  inverse  mapping  of  't(w),  then  we  can  write  the  Laurent  expansion  of 
(<M*))"  as 

(16)  (4>(*))n  =  7o*"  +  +  L  ^ .***• 

ksO  k=-o o 

The  principal  part  (or  polynomial  part)  of  ( <t>(z))n  is  defined  to  be  the  Faber  polynomial  of 
degree  n  for  the  set  fi.  It  turns  out  that,  on  assuming  some  additional  smoothness  on  the 
boundary  dfi  of  fi,  these  Faber  polynomials  can  be  used  to  generate  a  SIM  for  the  region  fi 
which  is  asymptotically  optimal  (i.e.,  an  AOSIM).  This  is  given  in  Theorem  22  of  [5], 

In  our  opinion,  these  results  from  complex  function  theory  form  the  foundation  of  a  wide- 
reaching  theory  of  iterative  methods. 


